Consistency in Formulation of Spin Current and Torque Associated with a Variance of 
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Stimulated generally by recent interest in the novel spin Hall effect, the nonrelativistic quantum 
mechanical conserved currents, taken into account of spin-orbit coupling, are rigorously formulated 
based on the symmetries of system and Noether' theorem. The quantum mechanical force on the spin 
as well as the torque associated with the variance of angular momentum are obtained. Consequently, 
the kinetic interpretation of the variances of spin and orbit angular momentum currents implies a 
torque on the "electric dipole" associated with the moving spin. The bearing of the force and the 
torque on the properties of spin current in a two-dimensional electron gas with the Rashba spin-orbit 
interaction is discussed. 
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PACS numbers: 72.10.Bg, 71.70.Ej, 72.25.-b 

The theoretical prediction of spin Hall effect (SHE)Q, 
H 01 has drawn considerable interest recently. The 
spin Hall current comes forward from the assembly effect 
that the carriers with opposite spins move in the oppo- 
site directions perpendicular to the electric field. Equiv- 
alently, a spin imbalance between its opposite flows is 
electrically generated via the spin-orbit (SO) coupling in 
the confined semiconductor. The manipulation of trans- 
verse force on the moving carriers with different spin 
orientations, i.e., pushes spin-up and -down carriers in 
the opposite directions, is understood as a dynamical is- 
sue of spintronics0, 0, Q. By using the optical meth- 
ods, the experimental detections of spin accumulation in 
connection with SHE have recently been reported 0, |t|. 
The effect of impurity scattering on SHE has also been 
investigated [13, 111 E3- 

However, a more fundamental question is emerged in 
the conventional understanding of SHE. The validity of 
the description for various physical properties of SHE de- 
pends on the correctness of definition for the spin current 
operator^, Q, 0|. Since the transnational and rota- 
tional motions are associated to the force and torque, 
without a rigorous definition of spin current, in some 
sense, the discussion about the convenient power source 
for the spin current is also somewhat arbitrary. Cur- 
rently, the spin current operator is practically defined in 
analogy with that of charge current as a tensor operator 
of symmetric multiplying of Pauli matrix and velocity, 
jf = (h/2){a a , vj [j, [1, [I . N evertheless it does not sat- 
isfy a continuity eauatioii|lllllll^. dpi /dt+ V- jf ^ 0, 
where = (h/2)^a a ip is a spin density. Unlike the 
charge current, the conserved spin current can not be ob- 
tained within the context of Schrodinger equation with 
the addition of the SO interaction. 

Actually, spin motion is intimately related to the or- 
bital angular momentum of spin carriers. If we identify 
spin as the internal angular momentum of carriers, the 
orbit angular momentum can then be regarded as a dy- 



namical property of carrier motion. The transmutation 
of spin orientation is compensated by that of orbit angu- 
lar momentum during the carriers moving in the system. 
Although the spin current is not conserved, the total an- 
gular momentum consisted of both spin and orbit angu- 
lar momenta should be conserved if there is no additional 
external force and torque. 

In this paper, we inherit the Noether's idea0|, and 
provide a systematical formulation of the conservation 
laws for the translational and rotational motions of spin 
carriers based upon the generic symmetry principles. 
Considering a system with Dirac fermions coupled to an 
electromagnetic field, we formulated the practical forms 
of spin and orbit angular currents with the symmetries 
being maintained in the sense of the validity of Noether's 
theorem. The nonrelativistic approximation (NRA) of 
conserved currents is derived with the conservation laws 
being guaranteed. Under this formalism, the torque on 
the angular momentum and the force on spin carriers can 
be identified rigorously. As a result we gain an insight 
of torque on the transmutation of spin orientation and 
orbit angular momentum, and ascertain the role of SO 
coupling in the conservation laws. 

Consider first the relativistically covariant Lagrangian 
of a closed system involved Dirac fields 'J', ^, and 
electromagnetic field A^. The corresponding La- 
grangian is given by £ = — (he/ 2)^^^ V ^ — mc 2 ^^ — 
(l/4)^ r At!y ^ r Aty 0, 01 , where the four- vector notations 
have been used (// = 1,2,3,4): x M = (x, ict), A^ 

I) -ia a \ ( I 



(A, itp), and 7 M 











(a = x, y, z) are Pauli matrices; 
(Vlf)g, V^d^~ i{e/hc)A„ 
operator and J-^ = d il A v 
strength tensor. From the Euler-Lagrangian partial dif- 
ferential equation, we can obtained the Dirac equation 
and the Maxwell's equation, 'y^V^ + (mc/h)^> — and 



/ p,h = .m„<i - 

is a covariant derivation 
d v A^ is the Maxwell field- 
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d v Fp.v = ie$l"ly$l, respectively. 

According to Noether's theorem, any continuous one- 
parameter set of invariance of Lagrangian is associ- 
ated with a local conserved current. For a closed sys- 
tem involved Dirac field and the electromagnetic field, 
the symmetries of system yield the conservation laws 
: (1) The U(l) symmetry leads to the charge conser- 
vation. Under the local U(l) gauge transformations, 

* — > exp(-ieA/ftc)\P, * — > * exp(ieA/Hc), and -> 
A fi — dfj,A , Noether's theorem gives rise to the conti- 
nuity equation d^J^ = with the charge current den- 
sity J M = iec^-f^- (2) Spatial translational isotropy 
begets the conservation of linear momentum while a time- 
translation symmetry gives rise to the energy conserva- 
tion. The continuity equation d^Q^ = is obtained 
from the invariance of Lagrangian under the time-space 
translation transformations Sx^ = ^ M „a„, <5\P = , 8^ = 

0, and 8A„ = 0. 9^ = {hc/2)^ < V ^ + T w T vp - 
\ T^i v/J 7 ^ v'S^v + dpiT^pAv) is the energy-momentum 
current density tensor; (3) Rotational isotropy gives birth 
to the conservation of angular momentum. The La- 
grangian is invariant under the transformations x'^ = 

{^tiu + e^Xu, tf'(a:') = [I-i{\/A)e llv (j llv \^{x),^ l {x') = 
ty(x)[I + i(l/4)e M „<7 M „], and A^x') = (5^ 

V ~1~ ^- 111*) A V {X) i 

with Gp V = i(l/2)[7 AI) 7 ! y]. The total angular momentum 
current density tensor J\ iliv = { x^Qxi, + J-\^A U — (p, <-> 
v)} — i(hc/4)^{a^, 7a}* can be defined to satisfy the 
continuity equation d\J\ tllv = 0. 

The NRA of Dirac equation can be car- 
ried out by separating the Dirac wavefunc- 
tion into large and small components[17| , i.e., 

* = (i/ja, ^b) T exp[— i(mc 2 + e)t/h]. Following the 
procedure of NRA given in Ref. [17], the normalized 
two-component wavefunction can be written in the form 
of tp = [1 — (ti 2 /8m 2 c 2 )a ■ Da ■ T)]ipA with a notation 
D = V - i(e/hc)A. The NRA Hamiltonian reads fiftKh} 
H = -(h 2 /2m)T> 2 + eip + i(eh 2 / 'Am 2 c 2 )a • (E x D) - 
(eh/2mc) a B (efi 2 /8m 2 c 2 )V • E - i(eh 2 /8m 2 c 2 )a ■ 
(V x E), where E = -Vip and B = V x A are the 
electric field and the magnetic field, respectively. 

Alternative to obtain the continuity equation from 
Hamiltonian we engage the same procedure of NRA to 
the relativistic continuity equations. The correspond- 
ing NRA of currents can be obtained from the rela- 
tivistic continuity equations by the power expansion in 
v 2 j <? with the same procedure as obtaining the Hamil- 
tonian with a SO coupling term. The symmetries of sys- 
tem are not broken by this power expansion. It is in- 
structive to demonstrate it in the NRA of charge cur- 
rent J M = iec^-f^fy , which satisfies d^J^ = in the 
relativistic quantum mechanics. In the nonrelativisitic 
limit to the order of v 2 /c 2 , the continuity equation be- 
comes p c + V • J c = 0. The charge density p c = 
eipty + (h 2 e/8m 2 c 2 )V 2 {^ip) - i{h 2 e/8m 2 c 2 )V ■ (^a x 
Di/>) and the current density J = —i(he/2m)i/j' Dip + 



(e 2 ft/4m 2 c 2 )[E+(2mc 2 /e)V] x^atp have been expressed 
in terms of the basic local observables. Look at the charge 
density first. The second term concerns the variance of 
localized density of carriers. The last term is in connec- 
tion with the accumulation of localized spin due to the 
spin flip. These two terms give a polarization vector 

P = -^[V(^)-i^xD^]. (1) 

This shows that the nonuniform and the SO interaction in 
the system contribute an effective electric dipole. In the 
charge current the first two terms are recognized as iden- 
tical to ep c v with the Hamiltonian given above, where v 
is the velocity vector. The last term is a "spin magnetiza- 
tion current" cV x m arising from an "intrinsic magnetic 
moment", where m = (eh/2mc)(ijj^aip). The supplement 
of a magnetization term in the charge current has been 
established for a long time[lf|. The Dirac equation im- 
plies the conservation law for the current d^J^ = as 
well as the decomposition = G + d^m^, where 
J^ G is Gordon current and m^ u — i(eft/4mc)$[7 /J ,7„]$ 
is interpreted as the magnetic moment density due to the 
carrier spin|l9|. 

The above analytic approach can be applied to the 
momentum current tensor 9 MJ/ . From d^Q^a = we 
have the continuity equation p M + V • $ M = 0, where the 
momentum density is given by 

+<# - 5%(Vp C ) x B + -E x B (2) 
8m z c 6 c 

with p^ = -i(h 2 e/8m 2 c 3 )^(<TB ■ D - a ■ BD)tp and 
= (fr/l6m 2 c 2 )V • OV x DDV) - (ft 2 e/8m 2 c 3 )V x 
(tp^tpB) + (l/c)V • (EA). I M is the nonrelativistic limit 
of momentum current which takes a form as 

4m 8m z c z 
+ O f CT x D V>)E - VO + V0E + (VE)^V] 

1 he 

-(EE + BB) + -(E 2 + B 2 B ■ V>W)I 

2 mc 

he h 2 <— > 

H Bip'aib + V x (BA - i — ib>a D V) (3) 

2mc 4m 

with the notation tp^ DpD a ip — ijA DpDaip + 
{D* p D* a ^ - (D* a ^)(D^) - (D^)(D a if>). 

The NRA of conservation of total angular momentum 
current tensor J\ !fiu can be obtained in the same way. 
The local angular momentum density is found as 

P ^ = X X/) M +p s + -ExA+-^/B, (4) 

where x = x — xo with a constant vector xo. The first 
term is the conventional angular momentum with the mo- 
mentum density, while the last two terms in Eq. (4) 
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come from electromagnetic field. The second term is the 
intrinsic angular momentum density, associated with the 
local spin, given by p s = (h/2)[l-(h 2 /8m 2 c 2 )V 2 ]^a->Jj + 
Pso + Pg> where p% Q = (fi 3 /8m 2 c 2 )[2ipia'D 2 ip + (Dt0t) • 
cr(D?/>) + (D^^cr • {Dip)) is the contribution of SO cou- 
pling induced by the transmutation and relates to the 
surface terms V • (ip^a D -0) and V x (ip^ D ip). The local 
angular momentum density satisfies the continuity equa- 
tion p A + V • JJ" 4 = with the total angular momentum 
current density tensor 



x x [J 



M 



i V x O f crD ip)] 

4m 



h 2 + <-> 
i ip' D aip 

4m 



+i— Vx WVD^i xxl + fB- A)I-AB, (5) 
4m 

where I is unit tensor. Notice that the survival of the 
term adhered the constant vector Xo in a conserved angu- 
lar momentum of system is not a trivial matter. The con- 
tinuity equation of angular momentum is satisfied with 
the symmetry of the invariance under spatial translation 
being guaranteed. Let us emphasize that the present 
form of total angular current is independent of the choice 
of coordinate origin in the system. 

So far we have got the NRA expression of the conserved 
angular momentum current for a closed system including 
the fcrmion and the electromagnetic field. In general, 
the NRA conserved currents consist of three parts: or- 
bital, spin, and electromagnetic. We next give a decision 
in favour of the angular momentum current of fermionic 
system by taking the electromagnetic field as an external 
field. Because the conservation laws are equations of mo- 
tion for the observables, the force and the torque on the 
spin carriers, coupled to the external field, can be deter- 
mined by the continuity equations. It is noted V • [(EE + 
BB)-(l/2)(£ 2 + B 2 )I + Vx(BA)] = E/ 7 + (l/c)(J c x 
B), which corresponds to Lorentz force, in the continu- 
ity equation of momentum flux when the Maxwell equa- 
tion is applied to the time-independent situation. While 
the terms, in Eq.(3), -(h 2 e/8m 2 c 2 )[i(ip^a x D-0)E - 
VO f V0E + (VE)^j + (he/2mc) [BipUip - (B • ip^aip)!] 
can be regarded as the source of the force associated with 
the carrier's spin0, [^J. The continuity equation of mo- 
mentum current becomes p M + V • I M = F, where 
the primes are meant to indicate the fermionic partition 
of density and current. The momentum density vector 
and the momentum current density tensor are given by 
p M ' = -i{h/2)ip^ip and J M ' = -(h 2 /4m)ip^b~Dip - 
i(h 2 e/8m 2 c 2 )V x (tp^a'Dip) -i(h 2 /4m)V x (ip^a D ip), re- 
spectively. The total force is found 



The first three terms are the Lorentz force and the 
" Stern- Gerlach force" which relates to a magnetic field 
gradient VB. The forth term is new which has not been 
discussed previously. It relates to the fluctuation of local- 
ized spin current and produces a torque on the angular 
momentum. The continuity equation of angular momen- 
tum has the general form p Al + V • J^' = M, where 
J A ' is the total angular momentum current of electron 
which can be written in the form of a sum of three terms 
: + z(ft 2 /4m)V x (V^cr D ip x x). The orbit 



I Al = J 



a 



angular momentum current is given by 



= xx $ M/ -i V x (ip^aTtip)] 

4m 



(7) 



and the spin current is defined by 



J 0a 



h 2 <— > 

4- 1p^U a Df}1p 

4m 



h 2 e 



E~^ip. (8) 



The term on the right hand side of continuity equation 
is the torque M which relates to the total angular mo- 
mentum, not spin partition itself. We find 



M = xxFH p* xB 

mc 



8m 2 c 2 



O f cr x D?/>) x E. (9) 



There are several origins in the torque M. Besides the 
conventional torque caused by the force F, the second 
term is the torque in the fundamental equation of spin 
resonance theory and the third term is a torque in- 
duced by the spin current. Because we are ultimately 
interested in relating the induced torque to the effect of 
SO coupling, a fluctuation in the spin current AJ 5 = 
i(h 2 / Am)%p^ {a x D)^ is introduced or, equivalently 



Ac/ a — £af3yJ[3y 



(10) 



By using Eq. (8) the relation between A j£ and the spin 
current can be put in the form e Q , i a 7 AJ^ = — Jp a - 
This vector has its own significance 2 lj . It contributes 
a force if there is a spatial change of electric field and a 
torque on the spin when the carriers moves in an electric 
field. From Eqs. (6) and (9) an force AF and a torque 
AM induced by the SO coupling may be obtained, which 
can further be separated into the contributions from the 
accumulation of spin and the fluctuation of spin current, 
in a compact form 



AF 

AM 



e / 2c(VB)- (VE) 
2mc 2 V -2cBx -Ex 



P 

AJ 5 



(11) 



F = e(E 



8m 2 c 2 



V^E)^ + [J — cV x m] x B 



ftp Pp- p ^ ^. 



It is noted that AJ 5 is just the "electric dipole" < 1 >= 
— < (fi/2c)x(—ih'D/m) > with the intrinsic magnetic 
dipole moment fi = (fte/2mc)er. Physically, it is thinking 
to be a spin vector rotates around an axis with a velocity 
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—ih(D/m) to give rise a magnetic current circuit, equiv- 
alently, an " electric dipole" . The torque on the " electric 
dipole" is < 1 > xE as expected. 

Next, we apply the formulism to a two-dimensional 
electron gas (2DEG) with the Rashba SO interaction. 
This can be achieved by taking E = (0, 0,E Z ) and 
A = in our derivations. In the follow discussion the 
indices k and / are specialized to 2D indices while a still 
is indicated in 3D space. The Rashba Hamiltonian is 
then given by H R = (l/2m)p 2 - (X/h)(p x a y - p y a x ) 
with A = {h 2 e/ Am 2 c 2 )E z . The momentum current 
in Eq. (3) reduces to J^ 1 ' = -(h 2 /4m)^ V fe v\ip + 
ie kjz (X/2)[^ajViip - {h 2 /2m\)d 3 {^a z Viij)]. It is 
easy to prove that J^ 1 = and Fp = for the sta- 
tionary states. The total angular momentum current is 
given by 



jAi _ j 

J ka — J > 



Ji 



/,7 V k 



4m 



Vx^VD^x^. (12) 



It is shown that the orbit part of angular momentum 
current density, jg, = -(ti 2 /4m)e aj ix :1 (^ V; V fc -0) - 
i(X/2)S az x • (yffV^), is nonvanishing only for a — 
z. The corresponding spin current is given by Jf a = 
—i(H 2 1 '4m)i/)> cr a V k^ — £kaz(X/2)(ip'ip). It is shown that 
the divergence of spin current, in general, satisfies the 
equation p s + V • J 5 = T with the torque T x r y \ — 

i(\/2){^a z ^~ 



straightforward calculation shows VfcJ^ = —T z . This 
leads to VfeJ^.' = 0, while j£{ = Tj, i.e. the compo- 
nent of total angular momentum perpendicular to the x-y 
plane is conserved, which shows clearly how the transmu- 
tation between the spin and the orbit angular momenta 
occurs in the system with a SO coupling, while the in- 
plan components of spin current are not conserved by 
requiring a torque on it. It is worth to emphasize that F 
is nonvanishing in the presence of an in-plan electric field. 
The force gives rise to an extra conventional torque. In 
general, it can be shown a relation between the torques 
T and AM, acting on the spin itself and on the total 
angular momentum, respectively, 



AM : 



8m 2 c ; 



: ipi(a x E) x DV>- 



(13) 



We note that there is an essential difference between 
the force F defined from the continuity equation of mo- 
mentum current and the transverse force f on the spin 
for driving spin current. The latter is the covariant force 
on a spin defined from the equation of motion 



rav = — mV 



f. 



(14) 



where v is the velocity of a wave packet of spin- 
polarized electron. In a Rashba 2DEG it is found Vj = 
-(S 2 /2m) ¥ ;t < D^ _ (2X/h 2 )e zjk p s k + {l/m)c jkl V k pf . F 



is the corresponding tensor of velocity flux given by 
■Mi ( o\ lta\* ., 75 xhe secon d term in 



{2\/h 2 )e zlk Jf l . 



VO and T z = -i{\/2)(^a ■ V tp). A Q ) 



the velocity flux can be understood as a drift velocity be- 
cause of bringing the electric field to bear on the moving 
spin. The covariant force on a spin is found as 

fj = ^-^zjkTk + tjka^kTa- (15) 

Replacing T k = —(2m\/h 2 )J kz into the second term of 
fj, we find fj = (Am 2 A 2 / h 4 )e Z jkJi z - Because of the op- 
posite signs for the contrary components of spin-polarized 
spin current, the different spins experience the transverse 
forces in the opposite directions |y, |jfl. 

The main results of this paper may now be summa- 
rized as follows. The conservation law of NRA conserved 
angular momentum based on the Noether's theorem has 
been formulated systematically and rigorously. It is iden- 
tified with that the torques are generated on a moving 
spin through the "electric dipole" when the electromag- 
netic field is regarded as an external field. As an example 
of the application, we consider a Rashba 2DEG and show 
that a torque will develop on the spin. The transmuta- 
tion between the spin and the orbit angular momenta 
has been demonstrated explicitly. Furthermore, the co- 
variant force on the spin in the presence of SO coupling 
is clarified. 
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